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ABSTRACT 

We show that in realistic cases of accretion in active galactic nuclei or stellar-mass X-ray binaries, 
the Lense-Thirring effect breaks the central regions of tilted accretion discs around spinning black 
holes into a set of distinct planes with only tenuous flows connecting them. If the original misalignment 
of the outer disc to the spin axis of the hole is 45° < 8 < 135°, as in ~ 70% of randomly oriented 
accretion events, the continued precession of these discs sets up partially counter-rotating gas flows. 
This drives rapid infall as angular momentum is cancelled and gas attempts to circularize at smaller 
radii. Disc breaking close to the black hole leads to direct dynamical accretion, while breaking further 
out can drive gas down to scales where it can accrete rapidly. For smaller tilt angles breaking can still 
occur, and may lead to other observable phenomena such as QPOs. For such effects not to appear, 
the black hole spin must in practice be negligibly small, or be almost precisely aligned with the disc. 
Qualitatively similar results hold for any accretion disc subject to a forced differential precession, such 
as an external disc around a misaligned black hole binary. 

Subject headings: accretion, accretion disks — black hole physics — hydrodynamics — galaxies: active 
— stars: neutron 



1. INTRODUCTION 

Accretion discs arc common in astro physics on all 
scales from proto stars to AGN (sec e.g. iPringld 119811 : 
iFrank et al.1 120021 ) . Many treatments assume that the 
disc is aligned with the symmetry axis of the cen- 
tral object, although there is often no a priori rea- 
son for this. The first widely-studied case relaxing 
this restriction was the evolution of tilted discs around 
spinning black holes. Until recently the standard pic- 
ture of tilted disc evolution was that, in the regime 
where viscosity acts diffusively (technically, a > H/R), 
the inner disc would align or counteralign rapidly with 
the hole's spin, with a smooth warp to the still mis- 
aligned outer par ts. This is often called the Bardeen- 
Petterson effect (|Bardeen fc Pettersonl 119751 but note 
their equations do not cons e rve angular momentum; see 
IPapaloizou fe Pringld H98l IPringld [1991 IQgilvid [19991 
for detailed disc ussions of the correct equations; and 
King et aT1l2005l for t he possibility o f counteralignment) . 

In a recent paper, iNixon fe Kind (|2012f) showed that 
this evolution can be very different for large inclinations 
of the disc and spin, and/or low values of the dimension - 
lcss viscosity coefficient a (Shakur a fc Sunvaevl I1973D . 
Enforcing the connections imposed by conservation laws 
between the various components ( ' radial ' . and 'v ertical') 
of viscosity (|Qgilvidll99"9l) . INixon fc Kingj (|2012ft showed 
that the viscous torques in the disc may be unable to 
communicate the Lense-Thirring precession efficiently 
enough to produce a smooth warp. Instead one expects 
a sharp break in the disc plane between the aligned inner 
parts and the misaligned outer parts, connected only by 
tenuous rings of gas with inclinations changing rapidly 
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across the break. We shall see from our Eqn. [7] below, 
that even if the viscosity coefficients remained constant 
with warp amplitu de, the disc would stil l break for re- 
alistic parameters. lLodato fc Pried (|2010() show that an 
assumed break of thi s type remains stab le in 3D simu- 
lations of such discs. INixon et al.l (|2012fl show that in- 
clined, partially counter-rotating gas orbits within an 
accretion disc lead to cancellation of angular momen- 
tum and thus subsequent accretion. The mass flow rate 
through the disc can, for a time, be increased by large 
factors up to ~ 10 4 times that of the corresponding disc 
with zero inclination. 

These results suggest that sufficiently inclined discs 
might break, and that if the precession rate of the inner 
and outer disc differs enough, discs rotating in opposed 
senses might interact and produce dynamical mass infall. 
We consider these questions in this Letter. We ask 

1) for realistic parameters, can the Lense-Thirring 
torques exceed the local viscous torques in the disc? 

2) If they can, can we arrange the broken disc such that 
disc orbits counter-rotate? 

2. WHERE DOES THE DISC BREAK? 

Here we estimate analytically the radius at which the 
disc is likely to break. This will give us an idea of the 
parameters which lead to breaking, and whether it is a 
common event. To break the disc, the torque result- 
ing from the Lense-Thirring effect must overcome the 
local viscous torques, or equivalently, the orbits in the 
disc must precess faster than the viscosity can communi- 
cate the precession. To illustrate this point, let us imag- 
ine the two extremes. If the viscosity is dominant, the 
precession is communicated throughout the disc instan- 
taneously, and the whole disc precesses rigidly. At the 
other extreme where viscosity is negligible, orbits at dif- 
ferent radii precess at different rates and the disc must 
break into many distinct rings. The nonlinear connec- 
tion between effec tive viscosity coefficients, enforced by 
conservation laws (|Ogilvidll999D . tell us that once a disc 
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starts to break in this way, the v iscosity evolves so as to 
reinforce the tendency to break ([Nixon fe Kind[20ll) . 

To calculate the breaking radius for a given viscosity 
we can assume that the disc has no initial warp. Then 
we can consider the usual viscous torque, and to a good 
approximation neglect the more complicated physics of a 
warped disc. The azimuthal viscous force per unit area 
in the disc is proportional to the rate of shear Rdtl/dR 
(where R is the radial coordinate and fi is the disc an- 
gular velocity) and so can be written as 



The typical radius at which the disc breaks is given by 



.dfi 
~d~R' 



(1) 



where p, is the dynamical viscosity. The area of an in- 
terface in the disc is 2irRH, where H is the disc vertical 
thickness. So the viscous force acting in the azimuthal 
direction is given by 

dn dtt 
F v = 2nRHpR— = 2-kRvY.R—, (2) 
dR dR 

where we have substituted the dynamic viscosity for the 
kinematic viscosity (p = pv\ and used S = pH. The 
magnitude of the viscous torque G v is given by 



G v = \RxF u \ = 2nRvT,R 2 (-Q') 



(3) 



where the prime denotes the radial derivative. This is 
the usual viscous torque at the interface of two annuli 
in an accretion d isc (see iLvnden-Bell fe Pringld 119741 : 
iFrank et all 120021) . For near-Keplerian rotation (fi re 
GM/R 3 ) it becomes 



G„ = 3tw£ (GMR) 1/2 



(4) 



The Lense-Thirring precession induces a torque with 
magnitude 

G LT = 2nRH \fl p X L\ = 2TrRHn p Y,R 2 n |sin<9| , (5) 

where the Lense-Thirring frequency f2 p = 2GJ\ l /c 2 R 3 , 
the disc angular momentum density \L\ = £i? 2 £l and 9 
is the angle between the angular momentum of the black 
hole and the disc. We assume the disc breaks when the 
Lense-Thirring torque tears gas off the disc faster than 
viscosity can make it spiral inwards, which requires 



Glt ^ G v . 



(0) 



Making the standard assumption of a thin Keplerian 
a-disc (jShakura fe Sunvaevl 119731) and using J\ y — 
aGM 2 /c, the radius at which we expect the disc to break 
is given by 



R 



break 



< - 



J) *» 



(7) 



where R g = GM/c 2 is the gravitational radius of the 
black hole. 

Equation j7]) looks plausible; it makes sense that at 
large viscosity, low spin or small inclination angles the 
disc cannot be broken, i.e. i?broak < Rg- Conversely, 
for low viscosity, high spin and/or large inclinations we 
expect the disc to break into distinct rings at some radius 

-Rbreak > Rg- 



-Rbreak < 350i? g Isin 



flP/3fM 2/3 ^r 2/a ( H / R \ 

1 Vo.5/ Vo.i/ uo- 3 / 
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(8) 

where we have parameterized using quantities typical for 
AGN discs. This radius falls within typical discs, sug- 
gesting that inclined discs near spinning black holes are 
quite susceptible to breaking. We caution that at ex- 
treme parameters this simple argument may not suffice 
to predict the behaviour of the system, although we ex- 
pect the general behaviour to hold. In the next section we 
confirm the breaking of the disc with numerical simula- 
tions. These simulations are a preliminary investigation 
into this problem and we intend to follow up in more 
detail in future publications. 

3. COUNTER-ROTATION 

Let us consider an inclined disc that breaks under the 
action of a strong differential precession. Its inner and 
outer regions precess almost independently at different 
rates. The precession timescale is much shorter than the 
alignment timescale, which must wait for precession to 
induce dissipation. So both the inner and outer discs re- 
tain their inclinations to the black hole spin. The outer 
disc remains almost unmoved, while the inner disc (typi- 
cally a ring of radial width ~ H) precesses rapidly about 
the spin axis. If the angle 9 between the outer disc and 
the hole spin lies between ^ 45 — 135° the inner ring must 
form an angle 29 > 90° with respect to the outer disc af- 
ter half a precession period. The rotational velocities are 
now partially opposed. 

This configuration is similar to t hose adopted in the 
counter-rotating disc simulations in INixon et al.l (|2012[ ) 
and so must result in rapid accretion. We note the prob- 
ability that a randomly oriented accretion event lies in 
the critical range of inclinations is given by the fractional 
solid angle as cos(7r/4) i.e. re 70%. In other words, disc 
breaking and dynamical infall from counter-rotating ac- 
cretion flows must be common in active galactic nuclei. It 
is also common in stellar-mass X-ray binaries if the spin 
of the black hole (or neutron star) accretor is sufficiently 
misaligned with the binary plane. 

We report two simulations of an inclined disc 
around a spinning black hole using the Smoothed 
Parti cle Hydrodynamic s (SPH ) code PHANTOM (see 
e.g. | Price fe Federrathl 
INixon et all 120121: iNixor] 



201(1 ILodato fe Price! 120101 
20121 ). SPH perform s well in 



modelling warped discs (jLodato fe Pried I2010D finding 
excellen t agre ement with the analytical treatment of 
IQgilvid (|1999f) . This is to be expected as both treat- 
ments solve the Navier Stokes equations with an isotropic 
viscosity. The c onnections bet ween the viscosity coeffi- 
cients derived bv lOgilviel (119991 ) therefore naturally holds 
in our numerical treatment. In the simulations reported 
her e we implemented the Lense -Thirring effect, follow- 
ing INelson fc Papaloizoul (|2000t ). The simulations use 
a disc viscosity with Shakura & Sunyaev a ~ 0.1 (cf. 
ILodato fc~P rice 201(3), a disc angular semi-thickness of 
H/R ~ 0.01 and the black hole has a spin a = 1. Initially 
the disc has no warp and extends from an inner radius of 
50i? g to an outer radius of 250i? g , with a surface density 
profile S = So(i?/i?o) _p arid locally isothermal sound 
speed profile c s = c St o(R/Ro) ~ q where we have chosen 
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Fig. 1. — Full 3D surface rendering of the small inclination sim- 
ulation. The whole disc was initially inclined at 10° to the hole 
with no warp. This snapshot is after approximately 500 dynamical 
times at the inner edge of the disc (50i? g ). 
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Fig. 3. — The accretion rates with time for the small (black) 
and large (red) inclination simulations. The accretion rate is in 
arbitrary units. The time unit is the dynamical time at R = R g , 
i.e. the dynamical time at the inner edge of the discs is ~ 350. 
Note this is the mass flow rate through R = 50-Rg and not the final 
accretion rate on to the black hole. 




Fig. 2. — Full 3D surface rendering of the large inclination sim- 
ulation. The whole disc was initially inclined at 60° to the hole 
with no warp. This snapshot is after approximately 500 dynamical 
times at the inner edge of the disc (50i? g ). 



p = 3/2 and q = 3/4 to ac hieve a uniformly resolved 
disc (jLodato fc Prin glc 20071) . The disc is initially com- 
posed of 2 millio n particles, which for this setup gives 
(h) /H « 0.8 (cf. lLodato fe Prl^l20T(ih . The two simu- 
lations differ only by the relative inclination angle to the 
black hole. 

The simulation shown in Fig. [T] has an initial incli- 
nation of 10°, and thus i?break ~ 40i?g, i.e. at a radius 
inside the disc's inner boundary, and so we do not expect 
the disc to break. This agrees with the simulation, which 
shows the usual (Bardeen-Petterson) evolution with a 
smooth warp. In contrast, the simulation shown in Fig. [2] 
has an inclination of 60° and therefore i?broak = H0i? g , 
i.e. we expect the disc to break. The simulated disc 
does indeed break, producing multiple distinct rings of 
gas with large relative inclinations. This leads to phases 
of strong accretion when the rings are highly inclined, 
and quieter phases when they are nofl The accretion 
rates for the two cases are shown in Fig. G3 

We note that this evoluti on is more extreme t han the 
ID simulations reported in iNixon fc Kind (|2012f ). How- 
ever the numerical method used there assumed the gas 
always remained on circular orbits, evolving purely by 
viscous diffusion. This was appropriate for the problem 

4 Movies of the simulations in this paper are available at 
|http: / / www. astro. Ie.ac.uk/users/cjnl2/tearing. shtml| 



studied there, namely, whether such a relatively orderly 
disc could break at all, gi ven the nonline ar evolution of 
the viscosity predicted by lOgilvid (|1999f ) . The current 
paper studies a dynamical problem, where the inclina- 
tion of disc orbits can change so rapidly that viscous 
diffusion of gas in circular rin gs is no longer an adeq uate 
approximation. Interestingly, INixon fe King! (|2012| ) did 
remark that their simulated disc rings appeared to be 
trying to break into more than one distinct plane. 

4. DISCUSSION 

We first discuss the possible arguments against this 
picture. The main unknown in this work is the nature 
of the viscosity controlling angular momentum transport 
in the disc. In this paper we have assumed that this 
can be modelled as an isotropic Shakura & Sunyaev a 
viscosity. There is a strong basis for assuming that the 
radial transport of angular momentum (gove rned by the 
azimu thal viscosity) is limited to a ~ 0.1 (|King et al.l 
|2007| ) and we expect discs around black holes to be very 
thin away from th e immediate vicinity of the strongly 
accreting hole (e.g. IKing fc Pringleir2007l . and references 
therein). However, the nature of the viscosity is un- 
known. In reality t he local viscosity i s like ly to result 
from MHD effects (|Balbus fc Hawlevi fl99ll) . and may 
well be anisotropic. The azimuthal shear is likely to be 
secular, with gas parcels continually moving away from 
each ot her, whereas the vertical shear is probably oscil- 
latory ()Pringle!ll99l . This is suggestive of a favourable 
anisotropy where the vertical viscosity does not strongly 
oppose breaking, but the result is simply not known. The 
consistency requirements worked out by lOgilviel (|1999f ) 
for a locally isotropic viscosity show that in a strong 
warp, the viscosity trying to hold the disc together is 
likely to weaken. There appears no reason to suggest 
this differs for an anisotropic viscosity. 

Another possible complication is the thermal evolution 
of the gas. As the disc orbits do not all lie in the equato- 
rial plane, they must shock and heat up. In the simula- 
tions above we have assumed an isothermal equation of 
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state, so this extra heat is assumed to be radiated away 
instantly. This is reasonable, as the densities in black 
hole discs are high and cooling is likely to be efficient. 
However if the disc cannot cool on the local precession 
timcscale, it may heat up significantly so that Eq. [7] is no 
longer satisfied. In this case the disc may rapidly thicken, 
perhaps even becoming thermally unstable, and switch 
to a different mode of both accretion an d warp propa- 
gation. On the other hand, iNixon et al.l ([2012D found, 
in counter-rotating disc simulations using an adiabatic 
equation of state, that although the gas dynamics can be 
strongly modified by gas heating, the net result in terms 
of rapid accretion is similar. 

Otherwise there appears to be no obvious reason why 
this behaviour should be suppressed. We therefore ex- 
pect this to be generic to most cases of accretion on to 
black holes (or neutron stars, since the Lense-Thirring 
effect applies here too), and more generally for the evo- 
lution of gas discs in the presence of a strong precession. 

5. CONCLUSIONS 

We have shown that for realistic parameters a ran- 
domly oriented accretion event on to a spinning black 
hole is likely to form a disc which is susceptible to break- 
ing at a radius close to the hole. If the angle between the 
disc and the hole lies between ~ 45 — 135° the interac- 
tion of partially opposed gas motions is likely, and leads 
to cancellation of angular momentum and rapid infall. 

This quasi-dynamical form of accretion, which appears 
to be a generic consequence of randomly oriented accre- 
tion on to a black hole, significantly alters the standard 
picture of slow viscous accretion. There is nothing to 
prevent a succession of events where rings break off the 
inner disc edge and then precess independently of those 
inside or further out (see Fig. [5]). It is reasonable to think 
of this process as tearing up the disc in a chaotic way. 

We shall consider possible consequences of this picture 
in subsequent papers, but note several points here. 

1. Tearing the disc can lead to rapid gas infall, but 
the long-term rate of central accretion is ultimately con- 
trolled by the outer disc. 

2. In a stellar-mass binary system this means that 
tearing modulates a quasi-steady mass transfer rate. The 
modulation might have large amplitude if the disc/spin 
inclination is high. Even if the inclination is modest, 
there are likely to be observable effects which could in- 
clude quasiperiodic behaviour such as quasiperiodic os- 
cillations (QPOs). Still more effects can occur if the in- 
falling rings shadow the central X-ray source. 

3. By contrast, in active galactic nuclei, the viscous 
timescale of the outer disc may easily exceed a Hubble 
time, and no steady state is ever set up. Thus tearing of a 
significantly inclined AGN disc may promote significant 
accretion when the central black hole would otherwise 



not gain mass at all. 

4. The torques we have considered here arc all internal 
to the disc - black hole system. They cannot affect any 
conclusions concerning the global conservation of angular 
momentum or mass of this system. In particular, con- 
siderations of the long-term evolution of black hole spin 
through accretion remain unchanged, whether the ac- 
cretion is assumed to be cohe rent (e.g.. iVolonteri fe Rees 
20051 iBerti fe Volonter|2008D or chaotic (King fe Pringle 
2001120071: King et al.l (20081) . 

5. We can expect qualitatively (and sometimes quanti- 
tatively) similar effects for other cases where an accretion 
disc is subject to a f orced externa l differe ntial preces- 
sion. Most obviously, INixon et al.l (|2011bD have shown 
that the effective potential experienced by a disc accret- 
ing on to a misaligned binary, as is thought to occur 
when supermassive black holes are close to coalescence, is 
extremely similar to that caused by the Lense-Thirring 
effect. For initial inclinations near to co- or counter- 
rotation the disc respectively coaligns or counteraligns 
(for the subsequent evolution of prograde circu mbinary 
discs see ICuadra et al.ll2009t ILodato et al.ll2009L and fo r 
retrograde circumbinary discs see INixon et al.l I2011af) . 
However, disc tearing changes this picture and may well 
bring gas into the close vicinity of the holes on near- 
dynamic al timescales, and thus help with the last parscc 
problem (|Begelman et al.lll980[) as well as feeding prob- 
lems. 

6. In general any black hole has some spin, and in 
general any accretion disc plane may be inclined to this 
spin. To prevent any of the effects we have discussed here 
from appearing, the misalignment must satisfy 



i • at <r 3a H 
sm0\ < , 

1 1 ~ 4a R ' 



(9) 



which is extremely small for realistic parameters. For 
example, a moderately thick disc with H/R = 0.1, a = 
0.1 and a low spin a = 0.1 must be inclined by less than 
4° to avoid this process. We suggest that disc tearing, 
particularly in the inner disc, probably occurs in many 
if not most cases of black hole or non-magnetic neutron 
star accretion. 
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